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Abstract 

The calculation of a large family of four point functions of general BPS operators in J\f — 4 SYM 
is reduced to the evaluation of colour contractions. For ~ BPS operators Oa the explicit results 
at order g 4 for the function (OnO'iOnOz) are given up to n = 6. The OPE of the general result 
is performed up to the second order in the short distance expansion parameter. Two examples are 
given, in which the mixing of the operators in the intermediate channel can be resolved using four 
point functions computed by this method. 



1 Introduction 



The pedagogical usefulness of Af = 4 SYM has been widely evidenced in the recent years. This theory 
has been extensively studied both for its relevance to the AdS/CFT framework (for reviews see [1]) and 
for its purely QFT aspects, in which it can be viewed as the simplest interacting (super)conformal field 
theory (CFT), completely determined by the choice of the gauge group and known to be finite [2]. In the 
superconformal phase the spectrum of this theory is very rich, and it is arranged in representations of the 
supergroup SU(2, 2|4), that contains the four dimensional conformal group 50(4, 2) and the R-symmetry 
group SU(4:). The operators of the theory are thus organized in superconformal multiplets, in which there 
are both primary and descendant conformal fields. In each supermultiplet, the conformal primaries can 
be identified with their scaling dimension A = A^ ) + j(g 2 ), the spin (ji, J2) and the Dynkin labels \p, q 1 r] 
of the R-symmetry representation to which they belong. The unique quantum number that depends of 
the coupling constant g is the anomalous dimension r y{g 2 ) 1 that is the same for the whole supermultiplet. 
It is useful to think of Af = 4 SYM as a CFT with additional symmetries that further constrain the 
dynamics. From the Af = 4 SYM calculations we can infer general properties of ordinary CFT's. 

In a general four dimensional CFT, conformal invariance puts tight constraints on the correlation 
functions of the primary operators [3]. Two and three point functions are completely fixed, up to multi- 
plicative constants related to the trilinear couplings of the theory 1 . The freedom in the four point function 
are arbitrary functions of the conformal invariant cross ratios. In a CFT, these correlators are thus the 
first ones whose space-time dependence is not fixed by the symmetries. In a theory in which Operator 
Product Expansion (OPE) holds, one can relate the knowledge of these correlators to the anomalous 
dimensions and the OPE structure constants, related to the trilinear couplings of the operators in the 
theory via the two point function normalization. The presence of an operator in an OPE channel is 
dictated by the fusion rules, that determine the conformal families of operators which can appear in the 
product of the two fields. Since the OPE implements conformal invariance, if one can obtain all the 
OPE structure constants and the anomalous dimensions of the operators of the theory, in principle all 
the n-point functions can be calculated. 

For these reasons the four point functions of operators in Af = 4 SYM have been extensively studied 
[4, 5, 6, 7, 8, 22, 10, 11, 12]. Particular attention was given to correlators involving protected operators 
that obey shortening conditions, and belong to short SU(2, 2|4) representations [13, 14, 15]. Correlation 
functions of such operators have additional non renormalization properties [16, 17, 18], that can be used 
to simplify the calculations. 

These functions exhibit short distance logarithmic singularities, whose coefficients can be interpreted 
in terms of the anomalous dimensions and trilinear couplings of the operators in the theory. Although 
the operators in the external points are protected, to their four point correlators contribute also non 
protected fields, that appear in the intermediate OPE channel. Hence, the study of four point functions 
of operators belonging to short multiplets can give us a lot of far from trivial information, also about non 
protected operators. 

In this paper we will show how to reduce the calculation at order g 4 of a wide class of four point 
functions of protected operators in a particular flavour projection, to the computation of some SU (N) 

1 The ratio of the square of the three point function coefficients and the product of the ones relative to the two point 
functions does not depend on the chosen operator normalization and can be interpreted as the invariant trilinear coupling. 
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color contractions. The described method is particularly efficient for four point functions of ^ BPS 
operators, due to the non renormalization properties of their (extremal) three point functions. However, 
we will see that it can be applied also to other operators belonging to BPS multiplcts. 

We perform also the OPE analisys, up to the second order in the short distance expansion parameter, 
of this class of correlators. If one is able to find a sufficient number of correlators in which the OPE 
intermediate channel contains the same operators, it is possible to extract the anomalous dimensions 
(up to order g 4 ) of the operators, together with the tree level and one loop trilinear couplings of these 
operators with the external ones. 

The paper is organized as follows: in section 2 we present our method, and we discuss under which 
conditions it can be applied also to non \ BPS operators. In section 3 we perform an OPE expansion 
of the general result, and derive a system of equations for the 0(g 2 ) trilinear couplings and the 0(g A ) 
anomalous dimensions of the operators in that channel. Then, in section 4, two examples of such a 
computation are given, in which there are enough correlators to extract the information for the operators 
involved. In the last section we comment our result and we give some outlook on how to further pursue 
the work following this approach. Computational details and other results are presented in the appendix. 

2 Four point function of protected operators 

Considering four point functions of protected operators presents several advantages. In fact, provided the 
finitcness of the theory, these four point functions are finite, and there are no subtleties in the calculation 
due to the renormalization of the operators in the external points. The logarithmic singularities we obtain 
by performing the OPE expansion depend only upon the short distance expansion parameter. 

In this section we will show how to calculate exactly, up to order g 4 , by simply performing color 
contractions, all the four point function of scalar operators of the type 



for a particular choice of the flavour representatives, in the TV = 1 formulation (see appendix). This is 
a four point function of four ^ BPS operators, two of them having dimension 2, and the other two with 
equal dimension A = n. 

Actually, we can go further, this method can be used also for calculate four point functions with other 
protected operators, not necessarily \ BPS. 

2.1 General setting 

Let's first define how we choose the operators in the external points. As explained before, we consider 
four scalar \ BPS operators of integer dimension A = n, in the [0, n, 0] representation of the R-symmetry 
group SU(A). In particular we will use two A = 2 and two A = n operators. The A = 2 operator is the 
scalar operator in the 20' representation, belonging to the ultrashort supercurrent multiplet. Written in 
the N = 1 formalism (see appendix) , the two representatives we choose are 



I /r\ A=n / -/-<A=2 //-) A=n / -/- ) A=2 \ 



(1) 




(2) 



(3) 
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In general, there are several \ BPS scalar operators in the [0, n, 0] representation. For our calculation it 
is convenient to make the following choice of the representatives: 

ol h3 =xt bcl ... Cn _ 2 r 1 4tt 1 ■■■t c r\ (4) 

+ 22:3 r> + a A.b 4.C1 j.Cj,_2 

o^ R ' =x? bcl ... Cn _ 2 4 4 4 ■■■4 > 

where the two color tensors X L and X R describe the color contractions which define the operators. In 
order to describe i BPS operators, they have to be totally symmetric, but not necessarily traceless. 
We have all the elements to define the four point function under consideration. 

G n (x 1 ,x 2 ,x 3 ,x A ) = (0 1 L v ' 3 {x 1 )Cl 1 (x 2 )0^ 22 '' 3 {x 3 )C22{x 4 )) ■ (5) 

Actually the function G n has also indices L, R that we omit for simplicity. This function can be expanded 
in powers of the coupling constant g, giving 

G n (x 1 ,x 2 ,x 3 ,x 4 ) = G ( °\x 1 ,x 2 ,x 3 ,x 4 ) + g 2 G£\xi,x 2 ,x 3 ,x 4 ) + g 4 G {2) (x 1: x 2 , x 3 , x 4 ) H , (6) 

and can be computed order by order in perturbation theory. The computation passes through the 
calculation of several Feynman diagrams, that can be organized in different classes. This classification 
will help us in developing an order g computation. From this result we can extract a lot of information, 
by performing the OPE and analyzing the short distance logarithmic singularities which appear. We will 
see in section 3 that for our purposes the non singular part in the short distance expansion at order g 4 
can be safely neglected, without lose of information. 

2.2 Classification of diagrams 

Let us now analyze the different contributions in the calculation of the correlator (5) at a generic pertu- 
bative order. The tree-level diagram of this four point function is depicted in fig. 1. Due to our choice of 




Figure 1: Tree level diagram of G n 

flavour representatives, there are three bundles of lines: from x\ to x 2 , from x\ to x 3 and from x 4 to x 3 . 

The first type of quantum corrections that we analyze are the ones represented by the diagrams in 
figure 2. This means that we collect in this group all the interactions that do not involve the x 4 — x 3 
or the x\ — x 2 bundle, respectively. The core of the diagram, i.e. all the interaction lines and the other 
spectators that one can have, is represented here with the dark bubbles. These quantum corrections 
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Xi X4 X\ X4 




X 2 x 3 x 2 x 3 

Figure 2: Quantum correction reducible to three point functions 

vanish. As an example, let us consider the first class of diagrams. Due to the presence of the spectator 
lines from X4 to x 3 it is as if the operator in the point x 3 is 

Ol% = 5"»X« Ci ... Cn _ 2 4 Cl ---4 Cn - 2 - (?) 

This operator has bare dimension n and belongs to the [0, n, 0] SU(4) representation. So it is \ BPS. 
Hence the calculation of these diagrams is equivalent to the computation of the correlation functions 

(Ol,(x 1 )Ol 22 '' 3 (x 3 )C 22 (x i )) , (0 1 L 1 -' 3 (x 1 )Cl 1 (x 2 )0^ R ,(x 3 )) , (8) 

where 0^ R , is defined in (7) and the operator 0\, arises as a contraction of O^ 1 ' 3 with C[ x : 

0\ l ^X L aacx ... Cn _^ c i---^ n -\ ■ (9) 

The correlation functions indicated in (8) are (extremal) three point functions of \ BPS operators, and 
therefore they do not receive quantum corrections [17, 18]. 

Technically, to define the two classes of diagrams written in figure 2, we must add the diagram in fig. 
3, that for the same reasons as above can be related to the two point function 

(01,(^)0^,(^3)) , (10) 

that has no quantum corrections. 

X\ X4 




x 2 x 3 



Figure 3: Quantum correction reducible to two point functions 

Hence all the non vanishing quantum contributions to the four point function in (5) arise from the 
diagrams in figure 4. Actually the second diagram class, depicted as it is, contains all the others, but 
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here we will interpret it as " all the diagrams that remain" . At a generic perturbative order, this class 
can be very large, but we will see that one can handle it rather easily up to order g 4 . 

We will calculate first G„\ so that we develop the machinery, then the calculation of Gn^ (and 
obviously Gn^) is simpler. 

2.3 Calculation up to two loop 

Let's analyze the first class of diagrams of figure 4. The diagram is equivalent to the computation of the 
correlation function 

\a-2 Med m b l(xi)C\ 1 (x 2 )4 C 4 d (x 3 )C 2 2(x i )) , (11) 

where we encoded in the tensor 

Y a b\cd = Xabci---c n -2 -^cdci---c„-2 (^) 

all the information about the specific operators we use in the calculation. This tensor is obviously sym- 
metric in a <-> b and c <-> d. The factor /fe^-:; , i L_ 2 is due to the spectator chiral lines connecting the 
points x\ and x 3 . So, for this class of diagrams, the entire information about the space-time dependence 
of the perturbative corrections comes from 

H ab \ cd (x u x 2 ,x 3 ,x 4 ) = (^(siJClite)^^)^^)) • (13) 

Let us analyze now the second class. Here there are all the diagrams that do not belong in one of 
the other classes. The interaction bubble must connect all the bundles depicted in the tree level diagram 
(figure 1). This is not possible at order g 2 , in which the interactions can connect with each other only 
two chiral lines. Therefore, at first order in perturbation theory this class is empty. At two loop it is 
equivalent to the correlation function 



(n-2) 2 (n -3)! 

(4^)«-3 (x? 3 ) 

where 



^3 Zabc\ d ef L abc \ def (xi,X 2 , X 3 , X±) , (14) 



Z abc\def = Xabc Cl ... Cn _ 3 ^defa--- c„_ 3 ' 

L abc \ de f(x u x 2 ,x 3 ,x 4 ) = 

(r i ^t(xi)Cl 1 (x 2 )4Ut e 4 f (x3)C 22 (x i )} - 5 l[ 2 H ab \ de (x u x 2 ,x 3 ,x 4 ) . 




x 2 xa x 2 x 3 

Figure 4: Interesting quantum corrections (note the x 3 <-» X4 exchange in the first class) 
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The term subtracted in eq. (16) belongs to the first diagram class in figure 4. Note that the relation 



abc\def 



Y ab \ de holds. 



We have reduced the calculation of a generic four point function G n at order g 4 to the calculation 
of two particular functions 2 , H ab ^ cd and L abc \ de f . I n these functions all the information about the space- 
time dependence is encoded. The other {g independent) information is in the two color tensors Y ab \ cd and 
Z a bc\def- To calculate Gffl we have to analyze the different color structures of the two functions H ab \ cd 
and L abc \ de f. Contracting these structures with the color tensors Y ab i cd and Z abc \ de f will provide the full 
result. 



2.3.1 Color structures at order g 4 

The computation of the two functions H ab \ cd and L abc \ d tf described in (13, 16) requires the calculation 
of many diagrams, which are drawn in the appendix. The functions have free color indices, so each of the 
Feynman diagrams is multiplied by a precise color structure. This structure has to be contracted with 
the color tensors Y ab \ cd or Z abc \ de f, depending on the class to which the diagram belongs. 

In table 1 we list the color structures we find, and the diagrams that provide it. Symmetrization of 



Structure 


Diagrams 


rjab\cd r r 
■"1 JiacJibd 


B2 , B 3 , B4 , B 5 ,B 7 ,Bq, B w , Bn , B u , _B 14 , B w , B 2 q 


j jdb\ cd £ £ £ £ 

n 2 = Jikajjkbjjicjied 


Ai , A4 , Bi , Bq , B$ , B12 , C\ , C2 ■ • ■ Cq 


rrab\cd _ f f f f 
^3 — J itaj jkb J jicj ikd 


B15 , Biq , Bn , Bis , B21 


rrab\cd f f f r 
n 4 — J lacj jkb J ikd J ijl 


A 2 ,A 3 


j- abc\def » „ r r 
l^l JkaeJ jbf Jijcjikd 


A 5 


j abc\def r r r f 
*- 1 2 Jjaejkcfjikbjijd 


B22 , B23 , B24 , B25 , B 2 6 , B27 



Table 1: Color structures and two loop diagrams 



the color indices is not taken into account. This is guaranteed by the contraction with the tensors Y ab \ cd 
and Z abc \ de j. Taking this symmetrization into account we can relate these structures. In fact: 



^ab\cdH 3 



ab\cd 



Yr 



ab\cd 



H. 



ab\cd 



N 



~2 Hl 



ab\cd 



^ab\cd^4 



-rj<xb\cd _ ^ \r 



ab\cd H'l 



ab\cd 



(17) 



2.3.2 Two loop result 



Given these formulae, one can write down the complete result at order g 4 for the two functions. For the 
function H ab \ cd we have 



^jab\cd 



+9" 



= g 4 N 8(A 2 + A 3 ) - 16(2 B 2 + B 3 + B 4 - B 5 + 2 B 7 + B 9 + B w - Bn 

+ Bi 3 + B14 — Bi 5 — B w — Bn + B w — B w — B 20 + B 2 i) i/° b ' cd 
8(2 A x + A 4 ) - 32(Bi - B 6 + B 8 - B 12 + B 15 + B w + B u + B 18 - B 2 i) 



+ 32(2 Ci + 2 C 2 - 2 C 3 - 2 C 4 + C 5 - C 6 + 2 C 7 + 2 C s - C 9 ) 



h: 



ab\cd 



(18) 



2 These are not correlators of well defined conformal operators, thus we do not have to expect them to be conformal 
invariant. 



The function L abc \ de f is 



j^abc\def 



16 A f 



T abc\def 



(19) 



The structure that multiplies L 2 is absent, since the linear combination of diagrams B 22 to B 2 -j which 
multiplies it sums to zero (see table 1 and appendix). 

Contracting these structures with the tensors Y ab \ cd and Z abc \ de f (coming from the contraction of 
the color tensors X L and X R , see cqs. (12), (15)) we get the full result, which will depend upon the 
contractions 



Y ■ Hi = Y ab \ cd H^ cd , Y ■ H 2 — Y ah \ cd Hl 
Putting everything together we obtain (see appendix) 



ab\cd 



N 



(n-2)! 



(4^)»*M>i,)»-^i 2 .^l4 B(r ' 5)2 



+(Y -Hi) $ (2) 



1 s 



r r 



s' s . 



■ $(2)( r 



Z • Li = Z abc \ def L 1 

{Y ■ Hi) (r + s - 1) + 2 
(Y .#!) + („ -2> 



afcc|de/ 



(20) 



(y ■ gg) 

TV 
(Z-Li) 



where r and s are the conformally invariant cross ratios 



™2 , 7 , 2 ~,2 ™2 

_ x 12 x 34 _ x 14 x 23 



. 71 2 , 7 t 2 



. 71 2 , 7 , 2 



, (21) 



(22) 



2.3.3 Lower order computations 

We performed the computation at order g . One can make exactly the same considerations to calculate 
these functions at order g 2 . As stressed before, we can do this using simply the first class of diagrams 
in figure 4, since the second does not contribute. In fact only one Feynman diagram will appear in the 
computation, multiplied by the structure Hi. So 



G£ ) (xi,x 2 ,x 3 ,x 4 ) = - 



7-4— r. 7-s l n 9 2 i ! 7 rff (r,s) (Y-Hi) . 



(23) 



Hence to calculate completely the function we need only the color tensor Y ab \ cd . 
Similarly, at tree level we find 3 



2)! 



Xi,X 2 ,X 3 ,X4) = 



Y ab \cdd ab d cd 



(4tT 2 )™+ 2 (Xl 3 ) 



(24) 



2.4 Application to more general operators 

All these considerations were done for 5 BPS operators. To neglect the contributions coming from the 
diagrams written in figure 2, we used the strong non renormalization theorem for the extremal three 
point functions of ^ BPS operators [18], valid at all orders in perturbation theory. Hence, for | BPS 
operators we can apply the same reasoning at all orders in g. 

At order g 4 our results actually can be extended to a more general class of operators. Our considera- 
tions are based on the total symmetry of the color tensors X L ' R defining the operators in external lines. 

3 We can use instead of Y ab \ cd the tensor Z abc \ de f, because we know that 6 c f Z abc \ de j = Y ab \ de . Our notation is more 
general, because when n = 2 the tensor Z does not exist. 
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A careful inspection shows that all the results would be the same provided that X L ' R = X^^,, ^ y 
i.e. with the symmetry only in the first two and in the last n — 2 indices. Note that with our flavour 
choice this is always guaranteed. Thus one has only to check whether the operator has such a flavour 
representative. 

In order to apply our method, the contributions coming from the diagrams in figure 2, must have 
neither quantum corrections at order g 2 , nor logarithmic divergencies at order g . So, if the operators 
Ol,r have vanishing one and two loop anomalous dimensions, we only need to check the vanishing of the 
order g 2 contribution to the three point functions 

(O l ,Q 2 o>O b ) , (O L Q2 >O w ) . (25) 

Here Ol',r> , coming from the contraction of two scalar fields that belongs to the original Ol.r are always 
^ BPS operators due to our flavour choice (see eqs. (9), (7)). This extends our technique also to other 
protected operators. It has been demonstrated that the three point functions of two \ BPS with one 
general BPS operator do not receive order g 2 corrections [19]. Hence our result can be extended to all 
the BPS operator that can be written as indicated in (4). We will give in section 4 an explicit example 
involving a \ BPS operator. 



3 OPE analysis of the four point function 



One can use the four point functions calculated in such a way to extract information about the theory. 
As an illustration, here we perform the OPE of the function in the OPE s-channel, in order to study 
the different contributions coming from the conformal families of operators that appear in the expansion. 
For this purpose we can consider 



G n (x2,x 3 ,X4,) = lim x-1 71 G n (xi,x 2 ,x 3 ,X4,) , 



(26) 



and expand this function in the limit X4 



X3. In this limit the final answer will depend on the three 



variables 



'23' Uj 2ii u '34 



, simplifying the space time dependence. 



This expansion will have logaritmic singularities, written in terms of the expansion parameter X34. 
In the OPE interpretation, the coefficients of these singularities can be interpreted as coming from the 
operators passing in the intermediate OPE channel, in particular their trilinear couplings (up to order 
g 2 ) and their one and two loop anomalous dimensions. If we have enough four point functions, these 
quantities can be obtained solving a system coming from the coefficients of the expansions. 



The expansion of G 



(0) 



does not have logarithmic singularities. It is simply 
(n-2)! Y ablcd 6 ab S cd 



G<p{x 2 ,x 3 ,x 4 ) 



T 2\n+2 



(47T j ■ ^34 

We shall write down the order g 2 and g A contributions, namely the functions G 
involve the functions B and & 2 \ given in the appendix. We start with 



(1) and Gi 2) , 



(27) 



that 



n m, , (n-2)! Y ■ Hi j 



-4- 2 A 



+ 



2A 



16 8 ' 

y-3 A 



• x 3 4 



x 



23 



(^23 ' X34) 5 



^23 



4 2. 
9 + 3 A 



+ - , (28) 



1 '23 
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lo (x 2 ) 

where the dots stand for higher order contributions in the expansion parameter x 34 , and A — ^j^l ) ■ 
There are two different expansions nested within each other: the X34 expansion, from which we recognize 
the bare dimension and the Lorentz structure of the operators that appear (see [20]), and the A expan- 
sion, that allows us to extract different quantities, namely the trilinear couplings and/or the anomalous 
dimensions. 

For the two loop contribution we get 



G£Hx 2 ,x 3 ,x A ) 



(n-2)! 



1 



+ 



X23 ■ X34 



+ 



finite tcrm+(^6 N(Y ■ H{) + 2(Y ■ H 2 ) + 3(n - 2){Z ■ Li)J A+ 
(N(Y -HJ + ^Y. H 2 ) + \{n - 2)(Z ■ Ll j) A 2 
hnite tcrm+(-2iV(y • - 3(Y ■ H 2 ) - H( n - 2){Z • Li))a+ 
(~N(Y ■ H X ) -{Y-H 2 )-\{n- 2){Z ■ L^) A 
finite tcrm+i (l7 N(Y ■ H x ) + 41(Y ■ H 2 ) + 40(n - 2){Z ■ L x )) A+ 
(i N{Y . Hl ) + H(y . H 2 ) + H (n _ 2 )(Z • L0)A 
finite tcrm+(l9N(T • ffi) - 32(y • if 2 ) - 49(n - 2)(Z- Li))a+ 
^(y-ffi) - 12(y-ff 2 ) - ll(n-2)(Z-i 1 )^)A 



(^23 ' ^34) 2 



1 

+ 36 



L 34 
c 23 



(29) 



We do not write explicitly the finite terms because their OPE interpretation contains also the g 4 part 
of the trilinear couplings, that will introduce new unknowns in our system. In general, to extract this 
quantity we need also other equations coming from three loop calculations [21]. 

From the leading terms in the expansion parameter (28, 29), one can see that in this channel the 
lowest dimensional operators that contribute in a non trivial way 4 have bare dimension A^ ) = n and are 
Lorentz scalars. The lowest dimensional rank 2 tensors that appear have bare dimension A^ ' = n + 2, 



due to the 



(3:23-3:34) 



terms. We will see at the end of this section how it is possible to extract information 



also from these considerations. 



3.1 OPE interpretation: the lowest dimensional intermediate operators 

As an example, we focus our attention to the operators of bare dimension n that appear in the OPE 
s-channel of this correlators. Their contribution is associated to the leading terms in the expansion 
parameter X34 in (28), (29), from which we can infer that they are all Lorentz scalars. Hence, in the most 
general case, a set of Lorentz scalars operators Oi of A' ) = n can appear. The OPE structure constants 

4 Namely with some quantum corrections, either in their scaling dimension or trilinear couplings. 
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of the operators in the external points Chq and Cn q l are defined as 

where CL,Rq% are the trilinear couplings with the operators in the external points and Mj are the two 
point function normalization of the intermediate operators: 



gWg) gx-y fy^-^ 

x 2n-2 (^j. _ y^2 y2 y x 2 y2 



The coefficients that appear in the OPE analysis do not depend the normalization of the intermediate 
operators. We can choose Mi(g) — 1, so that all the g-dependence of the two point function (32) comes 
only from the anomalous dimension. For this particular normalization the OPE structure constants 
coincide with the trilinear couplings. If we want to recover the latter for different normalizations we need 
to take into account equation (30). 

From the expansions (27, 28, 29) and the OPE interpretation [20], putting 

A 4 ( 5 2 ) = n + . g 2 7 f ) +. 9 4 7 f ) + ... , (33) 

C L , Rqi {g) = C^ Rqi + g 2 C^ Rqi + ■■■ , (34) 

we have 

E4«=0, (35) 

i 

1 vr (0) r (0) >y {1) - 2 ( w ~ 2 )' ry ha ^ 

LQi RQil1 ~ (47T 2 )™+ 3 ^ ^ ' 

E + = -^=^(y • Hi) , (37) 

I E ca4W 1)a - (fe^i K ' + ^ y ■ ^ + V n ~ 2){z ' ' (38) 

1 V^[^(0) ^(1) , r W r (0) \Jl), r (0) r (0) (2) _ 
~2^\S LQi RQl LQi RQl ' 11 LQi^ RQi ii - 

(39) 

= (b^Ti ( 6 ^ ■ ^1) + 2(T ■ F>) + 3(n - 2)(Z ■ L l} ) . 

We spend a few words about equation (35). The rhs is equal to zero because in eq. (27) we do not 
have a 5 1 2 term. This is the unique equation in which will also contribute protected operators whose 

x 34 x 23 

trilinear couplings do not renormalize. These operators do not appear at all in the other equations. Their 
contributions can be easily derived from a simple tree level calculation. 



3.2 Information from the OPE 

The machinery we developed allows us to write a system in which the unknowns are the anomalous 
dimensions of the operators that appear in the OPE channels, and the trilinear couplings of these operators 
(for a given choice of the normalization of their two point function) with the operators in the external 
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points. If we have enough four point functions, this system can be solved. The number of necessary four 
point functions depends on how many operators with the same bare dimension appear in that channel. 
As we know, this is determined by the fusion rules, which must obey also the constraints imposed by 
group theory. The representation to which belong the operators in the intermediate channel must be in 
the product of the representation of the external operators that select the OPE channel. 

One can extract information also in the opposite direction: inspecting the OPE expansion of G n there 
are some operators whose presence in the OPE channel is not forbidden by group theory, but they do 
not appear beyond tree level, or do not appear at all. We saw that in the s-channel there are no scalar 
fields (in the interacting regime) up to bare dimension n — 2, and no rank 2 tensor fields up to bare 
dimension n. When there are only \ BPS operators in the external points, the former is a consequence of 
the general non renormalization theorem for extremal three point functions. In fact one can see that in 
the product of [0, 2, 0] and [0, n, 0] at bare dimension A' ) = n — 2 only the representation [0, n — 2, 0] is 
allowed, hence in the internal channel only i BPS operators will appear. It is worthwhile to ask what is 
the particular set of operators that appears in our OPE expansion, since we can infer useful information 
also from the fields that do not appear. 

Before generalizing the considerations of this kind, we have to remind that we made a precise choice of 
the flavour representatives of the operators in the external points. This can project out some operators, 
that can appear for different external representatives. For an illustration of this see for example [22]. 
General conclusions can be drawn with the help of this method only after having checked that the absence 
of some field is not due to the particular flavour choice. 

4 Examples 

In this section we will perform the described calculation for the four point functions of the type G3 and 
G4 (sec eq. (5)), and see how these results can be used to extract trilinear couplings and anomalous 
dimensions of the operators that appear in the OPE channels. In particular we will study also a case in 
which the operators in the external points are not all \ BPS. In the case of only ^ BPS operators we 
calculate G n also for n = 5, 6. These results are given in the appendix. 

4.1 A case without mixing: the operator S 6 

A very lucky case is the one in which there is only one operator with a given bare dimension. In this 
situation we need only one four point function to obtain the desired coefficients. 

Here we will analyze the OPE expansion of the function G3, from which we can extract information 
about the operator 




(40) 



1=1 



This is a scalar operator of bare dimension A^ ) = 3, belonging to the representation 6 of the flavour 



group. 
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4.1.1 The four point function 

Let us analyze the correlator 

G 3 (x!,X2, x 3 ,x 4 ) ee (tr(0i0i0 3 )(a;i) Cj x (x 2 ) tr(^0^J)(x 3 ) C 22 (a:4)) • (41) 

This is a correlation function of two operators Q 2 o' with two 5 BPS operators, having A = 3 in the 
representation [0, 3, 0] = 50 of the flavour group. We call these operators Q50. In the language of eq. (4) 
we have 

250 3 = -^dabcfiifyl 1 ^50 ' = ~^abc<i>\ 4>\ 4>\ , (42) 

where d a b c is the totally symmetric tensor coming from the anticommutator of two generators of the color 
group SU(N), according to 

{T a , T b } = S ab — + d abc T c . (43) 

In the OPE between the operators Q^,o{x 3 ) and Qw'{x 4 ) we can find also the operator In fact 
from group theory we have that 

20' x 50 = 6 + 50 + 64 + 196 + 300 + 384 . 

We know that the first non trivial contributions to the OPE expansions of the functions G n come from 
operator of bare dimension A^ = n, so the contribution coming from <S 6 is the leading one in the 
expansion in X34. There are no other unprotected operators of bare dimension 3 belonging to the allowed 
representations. The other operator with such characteristics is Q50, that does not contribute beyond 
tree level. Then G 3 meets all the necessary requirements in order to calculate the anomalous dimension 
and the trilincar couplings of Se- 

As we saw before, to calculate the function G 3 we need to construct the color tensors Y ab \ cd and 
Z a bc\def an d contract them with the three independent color structures given above, see eq. (20). For the 
given four point function we find 

Y a b\cd = j^d abe d cde Z abc \ de f = —d abc d de f (44) 

V ^ = Y-H ^-Vf-* (45) 

Y g _ *(* a -l)(* a -4) z r _ *(* a -l)(* a -4) m 

2 64 1 64 ' v ' 

so 

G < f\x ll x 2 ,x 3 ,x 4 ) = , (47) 

,(D, ^_ (jV 2 -l)(7V 2 -4) r 

16(4^)6 xl 3 xj 2 xj 4 



G^(x u x 2 ,x 3 ,x 4 )=- K - ^/ A l 2 , 6 2 4 4 B(r,s) , (48) 



G 3 2 \x 1 ,x 2 ,x 3 ,x i ) = [ ^ A >X 7 l- T - r - r \ '-B{r,sf [r + s] (49) 

- <I>'"' ! -■- ! - <T> 1Z i ■ i ! T /-«J»>-(,-..s) } , (50) 



where r and s are defined in (22). 
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4.1.2 Anomalous dimensions and trilinear couplings 

After performing the expansion of the computed functions and the OPE interpretation we obtain the 
one and two loop anomalous dimensions and the trilinear couplings of <S 6 with Q 50 and 0,20'- From the 
general expansions we can extract 

C (0) 1 (N 2 -l)(N 2 -4) (1) 1 



^ ™- ^ 3N2 (52) 

Here the operator <S 6 is normalized such that: 



1 (N 2 - 1)(N 2 - 4) 1 

S W' <53) 

where A = 3 + 3 2 7^ + ff 4 7s^ H • F° r simplicity here we understood the insertion of the wave function 

renormalization. The anomalous dimensions derived here match [23]. 

4.2 A case with mixing: correlators of type G 4 

Let us consider the functions G4. The representation [0, 4, 0] of SU(4) has dimension 105. There are two 
operators belonging to this representation with dimension 4, that we call T> 105 and Q105, respectively 
the double trace and the single trace operator: 

VHf = (S ab S cd + 6 ac 5 bd + S ad S bc ) ^fy c 3 4 , (54) 
QioeP = ( tr (T a T b T c T d ) + symmctrization) tffyt&t ■ (55) 

So there are three four point function in this class, that we call Gw, Gqq and Gt*q- 

At dimension A = 4 there is also another BPS operator for which this method can be applied. This 
is the \ BPS operator 2? 8 4- In fact, one can easily verify that 

2 

N- 

is a good representative of this operator. So we can use for it our method, and build the correlators 
G d v, G d Q and G dd , where 2?84 is represented with d. Also the protected operator T>20' satisfies the same 
properties as 2?84) but we cannot use for it our method because the correct definition of this operator is 
less trivial, and has also g - dependent corrections [21]. 

In this way we can calculate the relevant contributions to the six four point functions under consid- 
eration. So we get the numbers written in table 2. 

4.2.1 Scalar operators of A^ ' = 4 in the representation 20' 

These four point functions contain enough information to solve the system coming from the OPE expan- 
sion, in which the operators that pass are the three operators d written in [24]. In fact, group theory 
gives: 

105 x 20' = 20' + 105 + 175 + 336 + 729 + 735 
84 x 20' = 20' + 84 + "70^ + +175 + 512 +729, (57) 

35+35 45+45 256+256 



T>H 33 = ( 5 ab 6 cd - 6 ac 6 bd - -f ace f bde ) 4>l<t>\<l>%<t>i (56) 
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Function 


Y ■ <r 


y 


Gdd 


(N 4 - 1)(N 2 + 1) 


7V(iV 4 - 1) 


G-dq 


J^(N 4 -l)(2N 2 -3) 


\{N 2 - 1){2N 2 - 3) 


Gqq 




1 ( a r2 i\/a r4 /■* a r2 i i o\ 

§37 (Ar - l)(iV 4 - QN Z + 18) 


s < 

Wd 


(1\ — 1 ) [1\ — 4 J 


iv (I\ — 1 J (7V — 4 J 


C ,„ 
<-s<ffl 


(i V — IJ^JV — 4 J 


j\r/ i\r2 i \ / i\r2 a \ 




1 (AT 2 1VW 2 A\(0 AT 2 1} 




Function 






Gdv 


N 2 (N 2 - l)(N 2 + 6) 


\N 2 {N 2 - 1) 


G-dq 


|iV 3 (iV 2 - 1) 


i7V(iV 2 - l)(N 2 + 6) 


Gqq 


1 AT 2 ! AT 2 IV Ar2 i «\ 


1 ( AT 2 1\( AT^ 1 ]\t2 i qo'i 
— 32 V _ / \ _ loiV oD J 


Gdd 


i(Ar 2 -l)(A^ 2 -4)(2iV 2 -3) 


-3(iV 2 - 1)(A^ 2 -4) 


GdD 


Ar 2 (7V 2 -l)(Ar 2 -4) 





GdQ 


|A^(7V 2 - l)(iV 2 -4) 






Table 2: Color contractions for functions of type Ga 



We shall analyze only the leading term in the OPE expansion parameter, that can be interpreted as 
contributions coming from Lorentz scalars with bare dimension A^ ) = 4. In the representations written 
above, with these quantum numbers there are the two \ BPS operator in the 105 (that do not appear 
beyond tree level), and six operators belonging to the 20' representation, two of them, K^ QI , belonging 
to the Konishi supermultiplet, while the other four are superprimaries. These are the protected operator 
2?20', and the three superprimary operators Oi. We have also two operators in the 84 with A^ ) = 4, 
namely T>s4 and /Cs4, another superdescendant of the Konishi operator. These operators in principle can 
appear only in the function Gdd, but it is easy to check that /Cs4 does not contribute while 2?84 gives 
only tree level contributions. 

So, at quantum level we have to deal with the six operators in the 20' representation. Actually only 
the three Oi pass in the expansions of these functions. In fact, as we told before, the V 2 o' operator is 
protected and its trilinear couplings with the external operators do not renormalize at order g 2 . The two 
Konishi superdescendants do not appear in this channel, either at tree level or at leading logarithms, due 
to the fact that they are bilinear in fcrmions at leading order. One can verify that this is true also for 
the order g trilinear coupling 5 . 

From these four point functions we can solve the OPE equations (35-39) for the three operators Oi. 
The results of this analysis are reported in [20] . 

5 Comments and Outlook 

We developed a technique which allows to calculate (at order g 4 ) the set of four point functions of 
protected operators defined by equation (1) by simply performing some color contractions. Particular 

5 Being bilinear in fermions these operators have odd g expansion in the trilinear coupling with the operators in the 
external points. 
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flavour representatives of the operators are considered. From these correlators, via an OPE analisys, one 
can extract the anomalous dimensions and the trilincar couplings of the operators that appear in the 
intermediate channel. Our considerations can be extended to a more general class of protected operators, 
in particular to those belonging to general BPS multiplets. 

We presented two explicit examples in which this method was used to extract information about the 
operators in the intermediate channel. We calculated the anomalous dimensions of the Lorentz scalar 
of bare dimension A^ ) = 3 in the representation [0,1,0] = 6, and of the three operators (9, in the 
representation [0, 2, 0] = 20' of bare dimension A^ ) = 4, together with their trilinear couplings with the 
external operators. Moreover, for | BPS operators the four point functions (1) are calculated also for 
n = 5 and n = 6. 

One may ask how restrictive is our particular flavour projection. It would be interesting to study this 
especially in the case of only ^ BPS operators. In other terms one has to find how many independent 
space time structures at order g A has this class of correlation functions. A preliminary analysis suggests 
that there are at least two independent structures. Hence in order to take them into account we have to 
consider at least two suitable choices of the flavour representatives. 
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Appendix 



Notation and conventions 

The field content of A/=4 SYM [25] comprises a vector, A^, four Weyl spinors, %p A (^4=1,2,3,4), and six 
real scalars, ip l (i=l,2,. . . ,6), all in the adjoint representation of the gauge group SU(N). In the W=l 
formalism the fundamental fields can be arranged into a vector superfield, V, and three chiral supcrficlds, 
$ 7 (7=1,2,3). The six real scalars, (p l , are combined into three complex fields, namely 

¥ = -5= (/ + , 4 = ^ - , (58) 

that arc the lowest components of the superfields <I> 7 and respectively. Three of the Weyl fermions, 
ip 1 , are the spinors of the chiral multiplets. The fourth spinor, A = ip 4 , together with the vector, A^, form 
the vector multiplct. In this formulation only a SU(3) x U(l) subgroup of the original SZ7(4) R-symmetry 
group is manifest. <I> 7 and <3?| transform in the representations 3 and 3 of the global SU(3) "flavour" 
symmetry, while V is a singlet. We can choose the normalization of the axial U(l) R-symmetry such 
that the three complex scalars (p 1 have charge +1, the vector A^ is neutral, the gaugino A has charge 
+3/2 and the spinors of the chiral multiplets ip 1 have charge —1/2. The I flavour index can be raised 
and lowered with Su, so we make no difference between (p 1 and (pi. 

The generators of the gauge group are taken in the fundamental representation with Dynkin index 
d(N) = |. To perform the contractions we use also the fusion and the fission rules 

tr (T a A) tr (T a B) = ± (tr (AB) - Itr (A) tr {B)j , (59) 

tr (T a A T a B) = ± (tr (A) tr (B) - Itr (AB)) . (60) 

We shall use Fermi-Feynman gauge as it makes the correction to the propagators of the fundamental 
superfields vanish at order g 2 [2]. Actually a stronger result holds, namely the anomalous dimension of 
the fundamental fields vanish also at g 4 . The terms in the action relevant for the calculation of the Green 
function we are interested in are 



So = j dxj d A 6 ($f <f>i - V a UV a ) (61) 

Sec = !^|^ 3 ^( e ^r bc ^^ (62) 
S cv = jdxj d 4 9 (2igf lmn <f> l ? V m <S> In - 2g 2 f lmn f n qr ^V m V^ Ir ) (63) 

Sw = jdxj d 4 6 (|/ abc [D 2 D a V a ] V b (D a V c )^j (64) 

where f a bc are the structure constant of the gauge group SU (N) . Due to the fact that all superfields are 
massless their free propagators have an equally simple form both in momentum and in coordinate space 
and we choose to work in the latter which is more suitable for the study of conformal field theories. The 
propagators of the chiral and the vector field are 

(^(xM^J^AM - |J e «a-*fc-*a><* i , (65 ) 



(VixiAAmxjJM) = -§t 6 -^, «»<») 

07T X i} 
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where 



Xij — Xi Xj 



On — Oi — 6j 



» dxt 



Computation of the relevant integrals 

The whole result, in the limit x\ — > oo (sec (26)) can be expressed in terms of two functions 

dx 5 



X i5 X % X k5 ^jk 



X s u \ X jk ' X jk I 



dx^dxQ 



2 „2 „2 „2 J2 — „2 * I „2 '2 ' 



x i5 x j5 x i6 x fc6 x 56 x jfe 

where the two functions B and <I>( 2 ) can be written as [26, 6, 4] 



B(r, s) =— | log(r) log(s) 



-2Li 2 



log 



r + s — l — y/p 



2Li, 



1 + r — s + y/pJ yi — r + s + y/p; 

$(2)( r) s ) =_L J 6 ( LU(-pr) + LU(-pa) ) + 31og(-) (Li 3 (-pr) - Li 3 (-ps)) + 

+ 1 log 2 (£) (Li 2 (-pr) + Li 2 (-ps)) + J logV) log 2 (p S ) + 

71-2 , / m / n 71-2 , 2/ s \ 7tT 4 1 

+ y log(pr) log(ps) + — log , 



with 



p = 1 + r 2 + s 2 - 2r - 2s - 2rs , 



P = 



OO u 

\z\<l ^ z fc 



1 - r - s + yfp 



Lijv(z) = 



k=l 



k N ' 



Both B(r, s) and <I>( 2 ) (r, s) are symmetric in r and s. Moreover they have singularities in (r = 0, s 
(r = l,s = 0) and (r = oo, s = oo). We will write now the expansion of these functions near 
singularities, 



e 2 /2 1 



a; 2 V 9 3 l0g V x 2 



+ 



(e-x) 2 /8 4 



a; 4 V9 3 



e 2 e-x 1 



e 2 (e-x) 2 /24 12, / e 2 
I log 



V 25 5 

We do not need the expansion of B(r, s) near the point (r = oo, s = oo) because of the relation 



s(r , s) = i B (I,£) = I B (V 

r \r r J s \s s 
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For the function <J>( 2 )(r, s) we have 



-i^il21-131ogl ? |-::l (IB 



2 /,2\ / g 2 



1 (e-x) 2 /367 /e 2 \ 2 

n — -401og U + 11 log 2 



O i O I 9 




(75) 



Feynman Diagrams 

Here we list the Feynman diagrams used in our calculation [7]. The results indicated are all in the limit 
x\ — ► oo, after multiplication by a factor of x\ for all the diagrams except A 5 and B 2 2 - B 2 7, that are 
multiplicated by x\. We indicate here also equalities between the diagrams that arc valid only in this 
limit. The complete x\ dependence of a correlation function (with well defined conformal properties) 
can be recovered after the substitutions x\\ — > x\^x\^ x\± — > %24 x i3i x 34 ~ * x 34 x i2- These results are 
grouped depending of the number of vector lines in it. Diagrams of type A are without vector lines, the 
type B contain one vector line and the type C two. 



^ = 74^4 /(4;2 ' 3) ' M = W?WJ {Wh ^=(4^ 2 ' 3 ' 4 )] 2 ' 



2(4^ 2 ) 8 :r 2 4 /(3;2 ' 4) ' M (4tt 2 ) 9 x| 4 



tf, ;::i:2.-n. ^ 5 = __ _ /(3; 2, 4) , (76) 

B\ = Bg = B 2 a = B21 = , B 8 = B4, A 3 = A 2 , B7 = £>i2 = -B14 = -B16 = — 2^ ' 

B 2 i=B 4 -7(l + £ --)A 4 , C 6 = -i-A 4 , C 7 = -^4 1; Sio = Bii = -^A 2 . (77) 

4 r r 4 r 4 2 

Here there are other relations, valid in the limit described above, that make some diagrams sums to zero: 

B5 = B 2 o = B 2 , Bq = B13 = Bu = Big = B3 , Bis = B15 , B23 = B25 , (78) 

Ci = C3 , C5 = — Cg = C4 , Cs = — C2 , B27 = B 2 q ■ (79) 
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X4 X\ X4 X\ X4 




x 3 x 2 A 2 x 3 x 2 Aa X 3 Ai 

X4 



x 3 



Table 3: Diagrams without vector lines. 



X4 Xi X4 X\ X4 X\ X4 




x 3 x 2 C 2 X 3 X ? C 3 X 3 x 2 Q 4 X 3 

X4 Xi X4 X\ X4 Xi X4 




x 3 x 2 c& x 3 x 2 q 7 X 3 X 2 (Js X3 

X4 

x 3 

Table 4: Diagrams with two vector lines. 
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Four point functions of higher dimensional operators 

Here we indicate the correlators G n , n — 5,6 with only \ BPS operators in the external points. Their 
complete form can be reconstructed from equations (6, 20-24). The SU(4) representations under con- 
sideration are [0, 5, 0] = 196 and [0, 6, 0] = 336. These operators are normalized so that their two point 
function is 

<Otol (0)> = 5? • (80) 

We indicate schematically (following the notation of eq. (4) , with complete symmetrization of the color 
indices) the definition of the operators, and the normalization of their two point functions: 



n — 5 



a -+ U*\ M (iV 2 -l )(^ 2 -4)(iV 4 + 24) 
Qi96 = tr (0 ) M Qlge = 



8N 3 

2 /IVJU2 



,,2\ 15(7V 2 -l)(A^ 2 -4)(iV 2 + 5) 
D 196 = tr (0 3 ) tr (0 2 ) M Pl96 = — ^ L A £ ; ( 8 1) 



Function 


y -<5 2 


Y-Hx 


^^196^196 


(Af 2 -l)(Af 2 -4)(W 2 +5) 2 
4JV 


5(Af 2 -l)(Af 2 -4)(AT 2 +5) 
8 


^X>196 Ql96 


(iV 2 -l)(iV 2 -4)(jV 2 -2)(iV 2 +5) 


5(iV 2 -l)(iV 2 -4)(jV 2 - 2 ) 


8N 2 


16JV 


^Ql96 Ql96 


(JV 2 -l)(N 2 -4)(JV 2 -2) 2 


(JV 2 - 1) ( N 2 -4) ( JV 4 +24) 


24JV 3 


1927V 2 


Function 


y -ir 2 


Z-L x 


^£>196 X>196 


7N(N 2 -1)(N 2 -4)(N 2 + 17) 


N(N 2 -l)(N 2 -4)(N 2 -23) 


16 


16 


^^196 Ql96 


(iV 2 -l)(W 2 -4)(9W 2 + 10) 


(iV 2 -l)(7V 2 -4)(JV 2 -30) 


32 


96 


^Ql96 Ql96 


(Af 2 -l)(AT 2 -4)(W 4 +24Af 2 -24) 


(JV 2 -l)(JV 2 -4)(JV 4 -24A r2 +72) 


3847V 


1152AT 



Table 6: Color contractions for functions of type G5 



n = 6 

45(7V 2 - 1){N 8 + 6N 6 - 607V 4 + 600) 

87V 4 ' 
135(7V 2 - l)(jV 2 -4) 2 (iV 2 + 8) 

30(iV 2 - 1){N 6 + 177V 4 - 72N 2 + 162) 
TV 2 

occn/(\r4 i\/j\r2 . n\ . 





^tr(0 6 ) 


M 6 




^tr(0 3 ) 2 


M 33 


G 42 


^tr(0 4 )tr(0 2 ) 


M 42 


q222 


^tr(0 2 ) 3 


Mq222 



(82) 
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Function 

Go 6 o 33 

Gq6Q222 

Gqzzqzz 

Gq33Q42 
Gq33Q222 
Gq42Q42 
Gq42q222 
Gq222 0222 



Y-8 2 

3(A f2 -l)(4Ar 8 ~26jV 6 + 155iV 4 -450jV 2 +450) 
327V 3 

9(jV 2 -l)(JV 2 -4)(2jV 4 -15JV 2 +30) 
32JV 3 

(iV 2 -l)(2jV 8 +41jV 6 -211iV 4 +630JV 2 -810) 
167V 3 

15(JV 2 -l)(iV 2 +3)(JV 4 -3jV 2 +3) 



2N 2 




2N 


27(Af 2 -l)(JV 2 -4)(7V 2 -8) 


9(Af 2 - 


-l)(iV 2 -4) 2 (JV 2 +8) 


32JV 2 




1287V 


3(A r2 -l)(A r2 -4)(A r4 +9N 2 -54) 


9(JV 2 


-l)(iV 2 -4)(jV 2 -6) 



16JV 2 

27(JV 2 -l)(JV 2 -4)(Ar 2 +3) 
6N 

(jV 2 -l)(jV 8 + 72jV 6 +47jV 4 -666A f2 + 1458) 
24JV 2 



JV 



Y Jh 

3(jV 2 -l)(JV 8 +6iV 6 -60JV 4 +600) 
1281V 3 

9(jV 2 -l)(JV 2 -4)(3jV 4 -2QjV 2 +40) 
1287V 2 

3(A f2 -l)(2jV 6 -13A f4 +60JV 2 -90) 
15(JV 2 -l)(jV 4 -3JV 2 +3) 



16JV 

9(jV 2 -l)(iV 2 -4) 
2 

(A f2 -l)(A f6 + 17A f4 -72jV 2 + 162) 
87V 



3(iV 2 -l)(AT 2 +3)(2Ar 2 -3) 
12iV(iV 4 -l)(iV 2 +3) 



Function 



Y-Ho 



Z- L 1 



Goeo e 

Go 6 o 33 

GqBQ42 
Gq&Q222 
Gq33Q33 
Gq33Q42 

^o 33 o 222 

Gq42 042 
Gq42q222 
Gq222 Q222 



3(jV 2 -l)(jV 8 +47jV 6 -324jV 4 + 1200jV 2 -1200) 
2567V 2 

9(AT 2 -l)(iV 2 -4)(iV 4 -20) 



64iV 

5(AT 2 -l)(JV 6 +5JV 4 ~18iV 2 -54) 
16JV 

3(iV 2 -l)(7JV 4 +2jV 2 -15) 
2 

9(iV 2 -l)(./V 2 -4)(iV 4 + 19Ar 2 -80) 
256 

3(iV 2 -l)(A r2 -4)(2Af 2 -9) 
8 

9iV(iV 2 -l)(iV 2 -4) 

(jV 2 -l)(jV 6 +53A r4 -39A r2 + 189) 
12 

AT(Af 2 -l)(7W 4 +62W 2 -45) 
12JV 2 (JV 2 -l)(iV 2 +3)(iV 2 + ll) 



3(jV 2 -l)(A r8 -35JV 6 +2047V 4 -1200iV 2 +2400) 
1024Af 2 

9(jV 2 -l)(jV 2 -4)(jV 4 -2QjV 2 +80) 
512JV 

(7V 2 -l)(jV 6 -64iV 4 +270jV 2 -540) 
64JV 

3(jV 2 -l)(jV 2 + 10)(2jV 2 -3) 
8 

9(jV 2 -l)(jV 2 -4)(jV 4 -lljV 2 + 16) 
1024 



64 

9jV(iV 2 -l)(jV 2 -4) 
8 

(iV 2 -l)(./V 2 +3)(iV 4 -58iV 2 +36) 
96 

jV(jV 2 -l)(jV 4 +53A r2 -18) 
4 

30JV 2 (Af 2 -l)(Af 2 +3) 



Table 7: Color contractions for functions of type Gq 
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